Mesoscopic admittance of a double quantum dot 
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We calculate the mesoscopic admittance G(uj) of a double quantum dot (DQD), which can be 
measured directly using microwave techniques. This quantity reveals spectroscopic information on 
the DQD and is also directly sensitive to a Pauli spin blockade effect. We then discuss the problem 
of a DQD coupled to a high quality photonic resonator. When the photon correlation functions can 
be developed along a RPA-like scheme, the response of the resonator gives an access to G(uS). 
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The possibility to couple nanoconductors to capaci- 
tive gates has been instrumental for exploring electronic 
transport in these systems. Applying DC gate volt- 
ages allows one to tune the energies of localized elec- 
tronic orbitals to perform the transport spectroscopy of 
a nanoconductor and reach various conduction regimes. 
Gates can also be coupled to AC electric fields, to ob- 
tain e.g. photo-assisted tunneling or charge pumping 
Recently, the mesoscopic admittance G(uj) of a single 
quantum dot subject to an AC gate voltage has been 
investigated experimentally [2]. The low frequency limit 
G(uj — » 0) c± —iujCmeso can be interpreted in terms of 
a mesoscopic capacitance C meso determined by the cir- 
cuit geometric capacitances but also by the dot energy 
spectrum, which sets the ability of the dot to absorb elec- 
trons. This problem has been discussed theoretically in 
the regimes of weak [3|, |4j] and strong Coulomb interac- 
tions [H, In a more quantum view, gates can me- 
diate a coupling between the electrons of a nanocircuit 
and cavity photons. This is widely exploited in the con- 
text of Circuit-Quantum ElectroDynamics. Coupling su- 
perconducting qubits to a coplanar waveguide photonic 
resonator allows an efficient manipulation, coupling and 
readout of the qubits @, E|. In the dispersive regime 
where a qubit and a resonator are strongly detuned, the 
cavity photons experience a frequency shift which reveals 
the qubit state. This shift is sometimes discussed in 
terms of the qubit mesoscopic capacitance [9|. The reso- 
nant regime leads to vacuum Rabi oscillations in which 
the nanocircuit alternatively emits and reabsorbs a single 
photon [10[. 

Double quantum dots (DQDs) are mesoscopic circuits 
which can be made out of e.g. submicronic two di- 
mensional electron gas structures [12j, or top-gated car- 
bon nanotubes[20]. These devices can be used to elab- 
orate various types of qubits [Ill4l3|. and offer inter- 
esting possibilities in the context of Circuit-Quantum 
ElectroDynamics [3, [Hj. The behavior of a photonic 
resonator coupled to a DQD has been recently studied 
experimentally [l6| . However, on the theoretical side, this 
problem has aroused little attention. Besides, the AC 
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FIG. 1: (a) DQD circuit considered in this article (b) Con- 
figuration used for the measurement of the DQD mesoscopic 
admittance (c) Coupling scheme to a photonic resonator 



gate-biasing of DQDs has been studied in the context of 
spin and charge pumping (see e.g. [17] and Refs. therein) 
and photo-assisted DC transport (see e.g. [18]), but no 
theoretical study has been performed in the context of 
mesoscopic admittance measurements. 

In the first part of this paper, we calculate the meso- 
scopic admittance G{uo) of a DQD. We show that this 
quantity displays a very rich behavior. In particular, it 
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is directly sensitive to a Pauli spin-blockade effect 1£ 
A measurement of G(oS) seems an interesting way to per- 
form the spectroscopy of a DQD, in the context of e.g. 
a qubit use, which can forbid invasive DC probes [Hj. In 
the second part of this paper, we discuss the problem 
of a DQD weakly coupled to a high quality photonic res- 
onator. The resonator could offer an alternative to direct 
AC gate biasing for measuring G(uj). When the photon 
correlation functions can be developed along a RPA-like 
scheme, both the dispersive and resonant behaviors of 
the resonator can be predicted from G(lj). We briefly 
discuss the range of validity of the RPA scheme in the 
non-interacting limit. 

We first discuss the mesoscopic admittance measure- 
ment (Figs, la and b). We consider two single-orbital 
dots L and R with orbital energies £l and £r, coupled 
together through a spin-conserving tunnel barrier with a 
hoping constant t and a capacitance C m . We note dd a 
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FIG. 2: Response function 11(0) versus the DQD gate voltages 
in the non- interacting/interacting cases [panels (a)/(b)]. The 
white numbers indicate the DQD most stable charge states 
[3. We have used a L = -0.1, a R = -0.5 and B = 0. In 
panel (a) we have used ksT /t — 0.1. In panel (b) we have 
used Ec (R) = E c , Um = 0.7E C , t = 0.1E C , and k B T = 0.02E C . 

the annihilation operator associated to an electron with 
spin a e {t,i} on dot d e {L,R}, h da = c\ a c d(T , and 
nd = not + n>di- Dot d is connected through a tunnel 
contact to a grounded reservoir, and connected through 
a capacitance Cg [C% c ] to a DC [AC] bias generator with 
voltage Vg [V ac (t)]. The reservoir states are described by 
annihilation operators c d k a - The full hamiltonian of the 
circuit writes (up to a term proportionnal to the identity 
operator), H x = H DQD + Hi + H ac withj2l| 

H DQD = Y\ (e d - (j[gn B B/2})fi dc7 + V\ n d (h d - l)E d c 

z — * a , a f * a 

+ U m h L n R + t ^ (cj^So- + h.c.) , (1) 

^ = H d>fc>(T {[^daCdka + h.C.] + edk<jc{ ka C dk<7 ^ , 

Hac(Vac(t)) = J2 d ^dfldVacit) , (2) 

e L(R) = Ei {R) [l - 2n L g {R) - 2nf L) (C m /C* {L) )] + £ L rm, 
rig = CgV g d /e and the total capacitance of dot d |22| . 

For later use, we define tunnel rates T d — "Kv^\t d \ /h 
to the leads, with v$ the density of states per spin 
for reservoir d. We note AA(t) = (A — (A)o) with 
(A)o the average value of an operator A for V ac = 0. 
From the linear response theory, one finds An d (uj) = 
e(a L Xd,L( UJ )^ a RXd,R( UJ )) V ac(oj) with charge correlation 
functions Xd,d'(t) = -i0(t)([h d (t),h d >}} . The charge 
of the capacitor plates connected to V ac writes Q ac = 
— qlltilc — aR?iRe + 2\2V ac . Therefore, one obtains 

AQ ac (uj)/V ac (uj) = 2A 2 - (e 2 /h)U(u) = G(uj)/(-iu) 

(3) 

with n(cj) = J2 d d ,OLd®-d'Xd,d'{u) and G{uo) the admit- 
tance of the DQD. The term in A2 corresponds to the 



DQD response for totally closed quantum dots (i.e. t d = 
t = 0). In the low frequency limit, i.e. uj much smaller 
than the characteristic energies involved in the DQD 
dynamics (including r^m), we obtain C meso — 2X2 — 
(e 2 /h)U(0) G R. One can calculate n(0) from the defini- 
tion of Xd,d'(t). Alternatively, assuming that til and fiR 
have finite correlation times, i.e. lim^ +00 Xd,d f (t) = 0, 
one can use 

n (°) = h ^Zd4' a d a d'd (h d ) Q /de d i (4) 

with, assuming that the eigenstates of Hdqd (with 
energies Ei), are thermally populated [18|, 

(n d ) = Ei(^i\n d \i; i )exp(-l3E i )/j: i exp(-(3E i ) (5) 

We first discuss the non-interacting limit, using 
Ec iR \u m -> 0, C m = 0, B = 0, and T L{R) = T, 
which yields e d = —CgeVg/C^. In this case, II (o;) can 
be expressed exactly as TL(jjj) = IIi(o;) + 112(0;), with 

n i[2](^) = E S G{ + ,-} n 5, S [ S l(^). 

n s ,s>(oj) = — / deYf(e)g s ^(uj)/[(e-E s ) 2 +Y 2 ] (6) 
™ J-00 

and g 8i8 ,(u) = 8 , (e - E 8 ,)/[(e - E s ,) 2 - (uj + iT) 2 }. 
Here, s denotes the sign opposite to s. We use A S)S = 
(a L + + s(a L - a jR )cos[6>])/2, A S) 5 = -(a L - 
a R )sm[6}/2, = arctan[2t/(e L - e fl )], J5± = (e L + 
± A c )/2, and A c = a/^ - e^) 2 +4t 2 . In the limit 
T = and cj = 0, we obtain IIi(0) -2^X) S A 2 s ^ s with 
^ s = r/7r[^ s 2 + T 2 ] the DQD partial density of states 
(DOS) corresponding to state s dressed by the leads. 
This result is reminiscent from the non-interacting sin- 
gle quantum dot case 0, [23| where the dot DOS plays 
a crucial role. The term 112(0;) is more specific to the 
DQD case and is not simply related to v s . It is finite 
when aL 7^ (%r, i.e. when V ac induces different renor- 
malizations of the levels cl and e#. Processes which 
involve electronic transfers between the two dots thus 
contribute crucially to 112(0;). We now focus on the limit 
< T < fc^ < At low frequencies, we obtain from 
Eq.® 

IMO) = -pny^ A^cosh- 2 [^ s /2]/2 (7) 

Z J s 1 

n 2 (0) = -Ah(a L - a R )H 2 (f(E.) - f(E+))/A 3 c (8) 

with f(e) = 1/(1 + exp(/te)). These results can also be 
obtained from Eqs.(j4]) and (j5j). Hence, 11(0) does not 
depend anymore on T. Figure [2^ shows 11(0) versus the 
DQD gate voltages. The weak resonant line crossing the 
Vg = V g R = point is due to 1^(0) while the anticross- 
ing lines are due to IIi(0). For lj, 2t ^> T, we obtain 
n(o;) - n 2 (o;) with n 2 (o;) ~ 4h(a L - a R ) 2 t 2 (f(E_) - 
f (E+ ) ) / A c (o; 2 — A 2 ) . Dot /lead electron transfers are not 
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Along line c, 11(0) reaches a maximum which is y/2 
higher than along line <i, because lines c and d involve 
resonances between a different number of states. For 
uo ^> r L (^), using a master equation approach, we find 
n(cj) - -Il(0)A 2 (d) /(cj 2 - A 2 c(d) ) along line c(d). The 
finite frequency behavior of H(uj) will be discussed in a 
more complete way elsewhere. 

We now discuss the effect of a Zeeman field B on 11(0) 
(see Fig. 3). We use B > so that f spins have a lower 
energy. Lines of type a or b are shifted by B because 
they now correspond essentially to a transfer of t spins 
between the dots and leads. However, their height is 
almost not modified (except too close to triple points). 
For a magnetic field gfisB ~ t, 11(0) cancels in a region 
where line d was formerly extending (see Fig. 3b). This 
is because in this area, the state (t, t) becomes the most 
stable state, and therefore charge fluctuations between 
the two dots become impossible. This effect represents a 
near-equilibrium version of Pauli spin blockade [19| . As 
a result, line d is shifted to higher (lower) values of 
(rig), and it reaches a higher maximum which depends 

L( R) 

strongly on T. Indeed, we obtain for q\±bB <C E c y , U m 



FIG. 3: Effect of a Zeeman field on 11(0). We have used 
gfiBB — 0.2E C for panel a and the red full line in panel b, 
and g\iBB — for the blue dashed line in panel b. The other 
parameters are the same as in Fig.2.b 



relevant anymore because they are too slow, and II (a;) 
shows a resonant behavior due to the internal dynamics 
of the DQD. 

We now discuss the interacting case for < ^l(r) <^ 
k B T <C t <C Ec {R \u m . The DQD stability diagram cor- 
responds to the standard honeycomb pattern [21[. Figure 
EJd shows the variations of 11(0) with n^ R \ calculated 
from Eqs. ((4]) and ((5]) for B = 0. Different kinds of reso- 
nant lines occur in this graph. The first kind corresponds 
to electron transfers between the DQD and a lead, and 
has a width set by T. For instance, line a corresponds 
to transitions between states (0,0) and (0,<r), with a G 
{t, 1} [HI- This line can be approximated (away from 
triple points) as 11(0) ~ — to^/3/4cosh 2 [/3e^/2], which 
is reminiscent from Eq.©. Similarly, line b corresponds 
to 11(0) ~ -Ha 2 L /3/4cosh 2 [/3e L /2]. The second kind of 
resonances corresponds to electron transfers between the 
two dots, in the same fiR + til subspace. For instance, 
line c involves resonances between DQD states (<r, 0) and 
(0, a), with (7 G {t, I}- It can be approximated by 11(0) ~ 
—2h(aL — aR) 2 t 2 /Al, which recalls Eq.®. Line d cor- 
responds to a resonance between (0,ti), (t?i) an d (i,t)- 
It can be approximated by 11(0) ~ —4h(aL — aR) 2 t 2 / A? d 
with A d = ^(E 02 -E 11 ) 2 ^8t 2 and E 02 - E n = e R - 
cl + 2E^ — Um. The above expressions again do not in- 
volve the values of the tunnel rates due to ^l(r) <^ ^bT. 
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(9) 

with A = (En — E20 + Ad)/2. In contrast, line c is 
not affected by a magnetic field qiibB ~ t. Line c is 
affected by B once (t, t) becomes the most stable state 
near — rig — 0.5, which occurs only for higher values 
of magnetic field q\±bB ~ J7 m (not shown). 

To conclude this first part, mesoscopic admittance 
measurements appear as an interesting alternative to 
charge sensing [20, 25], for performing the spectroscopy 
of quasi-closed multi-quantum-dot systems. We have 
mainly discussed the ^l(r) <^ fc#T limit. The frontiers 
between the different (til^ur) domains can be seen in 
11(0). In the interacting case, the parity of the DQD to- 
tal occupation number can be determined directly from 
the difference of amplitude between lines of type c and 
d obtained at B — 0, or from the spin blockade effect 
obtained for B ^ 0. The DQD mesoscopic admittance 
also gives a direct access to information on the DQD 
spin state, since spin singlet and triplet states can be 
discriminated using spin blockade. At high frequencies 
uo ~ t ^> T, Ii(uo) shows resonances due to the internal dy- 
namics of the DQD. We have disregarded spin and orbital 
relaxation effects (with rates denoted T s j£), which can be 
due e.g. to magnetic impurities, spin-orbit coupling, or 



phonons. However, assuming Y 
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< k B T, the 



results presented here [Eqs. ([7]) to (|5J)] will not be affected 
for oj much smaller or much larger than T 
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an intermediary value of cj, the expression of U(uj) can 
involve explicitly and ^l(r)- 

We now consider an experiment where the DQD is con- 
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nected through C^ c and C^ c to an external (L r ,C r ) circuit 
which is a simple model for a photonic resonator [7j. The 
full circuit hamiltonian H2 includes terms in (C r /2)V 2 C + 
(l/2L r )4>^ c and XiV ac + \2Vac due t° the resonator and 
DQD respectively, with 4> ac the flux operator through 
the inductance L r and V ac the operator associated to 
V ac - We define the charge operator conjugated to 4> ac as 
Qac = V a c/C' r with C' r = C r + 2A2 and the photon anni- 
hilation operator a = —i/^2hZ r ^ ac + ^/ Z r /2hQ ac with 
Z r = ^L r /C' r . We assume that the resonator photons 
are coupled to an external photonic bath corresponding 
to the annihilation operator b [26]. We finally have 

H 2 = Hdqd + eVrms ^2 d a d h d (a + a 1 " ) + Ai (a + a 1 " ) 
+ hjj' r a^a + 7 fjLJpPbp + (rWa + r*a^b p ) 
+ (^ae iu;dt + ^e-^^a 1 ") + iJ z (10) 

with uj' r = l/y/L r C' r and V rrn8 = ^hbj' r /2C' r . 
The terms in ft account for an external driving of 
the resonator at frequency o;^/27r|27j. For simplic- 
ity, we study the response of the resonator through 
its mean voltage. The linear response theory gives 
AV ac (t) = Re[G a+at , a t(^)ft*e-^ £ ] with G Afi (t) = 
—i0(t)([A(t),B]) K= Q. We can relate Gat, a an d Gat, at to 
Xd,d'(t) = — i0(t) ([fid(t), fid']) by using an equations of 
motion approach, which takes into account the station- 
narity of Gat [a], at- We assume that the self-energy terms 
zCp \ T p\ 2 / (hLU^zhuj p -\-iO + ) due to the coupling to the outer 
photon bath write — ihA, with A > 0, to account simply 
for the finite quality factor of the resonator. We ob- 
tain the exact relation G a ^ a (uj) = Go + Go^ ms Il(uj)Go 

with G = (cj - uj' r + zA)"\ U(u) = Y^d4' a d 0L d'Xd,d'{ UJ ) 
and uurms = V rms e/h. Using an analogous expression for 
G a t, at and assuming A <C uj' r , one finds G a+ at, at — G a ,at- 
To find t he p oles of Ga,at, a self-consistent approach is 
necessary (ill [Isj. We postulate a RPA-like approxima- 
tion G a ,atH = G + G ^ ms n(o;)Ga,at(^), which yields 

G7 1 at H = G - 1 -o; r 2 ms nH (11) 

In the limit where hu f r and hu 2 ms ll(0) are both much 
smaller than the energy scales involved in the DQD 
dynamics, Eq. ([TT]) gives a dispersive shift of the pho- 
tonic resonance frequency, i.e. cj*°* c± uJ r + co^ ms II(0). 
This result can be recovered by considering a classical 
parallel (L r ,G r ) circuit in parallel with a capacitance 
2A 2 - (e 2 /h)U(0) following from Eq. Indeed, as- 

suming uj 2 ms ll(0) <C u)' r , we expect free oscillations with 
a frequency (L r [C' r - (e 2 /ft)n(0)]) -1 / 2 ~ cj*°*. For larger 
values of uo' r , in the general case, the response of the res- 
onator is not simply given by Tl{uj' r ) but by the func- 
tional form of H(<u) [and thus G(uo)). For instance, let 
us use the resonant form II (a;) c± Q/(uj 2 — A 2 ) obtained 
previously. One expects an anticrossing effect when the 



photonic resonator becomes resonant with the DQD. 
From Eq. ([Til), we indeed obtain = (A + u' r ) /2 ± 
v/A+(A-o;;) 2 /4 with A = (^o; 2 m J/(A + u' r ). 

In the non- interacting case, the RPA-like approxima- 
tion of Ga at can be justified by using a standard dia- 
grammatic perturbation theory in oll^r)- For each order 
in a L ( R ^ : the contribution to G a ,at corresponding to a 
series of "bubble" diagrams must be dominant. In prin- 
ciple, an estimation of diagrams at fourth order in a L ^ 
already provides a good indication on the validity of the 
RPA scheme [30|. From a dimensional analysis, the RPA- 
like development of Ga,at is valid at least in the regime 
T = with A,E±,huj' r ,huj t r ot - hw' r <C T. Considering 
the relevance of the results given by Eq. ([TT]h the RPA 
scheme is probably valid in a much wider range of pa- 
rameters. However, from the fourth order diagrams, it 
seems crucial to have u;£ ot — uo' r and A small, and V finite, 
this assertion b eing difficult to define quantitatively in 
the general case[3lj. 

As a conclusion for this second part, we have discussed 
the behavior of a high finesse photonic resonator coupled 
to a DQD. When photonic correlation functions can be 
developed along a RPA-like scheme, both the dispersive 
and resonant behaviors of the resonator reveal informa- 
tion on the DQD admittance. 

We acknowledge fruitful discussions with B. Dougot. 
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